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Abstract 

We study the discrete causal set geometry of a small causal diamond in a curved spacetime 
O using the average abundance (Ck) of A:-element chains or total orders in the underlying causal 

set C. We begin by obtaining the first order curvature corrections to the flat spacetime 
?-h . expression for (Ck) using Riemann normal coordinates. For fixed spacetime dimension this 

allows us to find a new expression for the discrete scalar curvature of C as well as the 
time-time component of its Ricci tensor in terms of the (Ck)- We also find a new dimension 
estimator for C which replaces the flat spacetime Myrheim- Meyer estimator in generic curved 
spacetimes. 

1 Introduction 

That there is a profound relationship between order and Lorentzian geometry has been evident 
ever since the work of Malament, Hawking and others [U [2] where they showed the existence of 
a bijection between the causal structure, itself a partially ordered set, and the conformal class 
of the spacetime metric. This is one of the main motivations for the causal set approach to 
quantum gravity, which assumes that the primitive structure underlying spacetime is a locally 
finite partially ordered set, or causal set [3]. Instead of considering the spacetime metric as the 
fundamental dynamical variable in causal set theory (CST) it is the causal structure that one 
wishes to "quantise" . However, to recover the full spacetime geometry from the causal structure, 
there must be a way to obtain the spacetime volume or equivalently, the conformal factor. This 
is achieved in CST via the condition of "local finiteness" which implies a fundamental spacetime 
discreteness: underlying every finite volume region of spacetime is a finite cardinality causal set. 
Thus, the continuum-discrete correspondence in CST is not exact but approximate, with the 
continuum being the approximation of the underlying causal set. 

In order to maintain the relationship between volume and cardinality in all coordinate sys- 
tems, a causal set C approximated by a spacetime (M, g) is obtained via the following random, 
Poisson discretisation of (M,g) [31 Q]. Given a fundamental scale p~ l (which could be the Planck 
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volume) the probability that a spacetime region of volume V contains N elements of C is given 
by the Poisson distribution 

Py(A0=exp-^^^, (1) 

for which 

(N) = pV, (2) 

thus establishing the required number to volume correspondence. To give credence to the exis- 
tence of a fundamental spacetime discreteness CST moreover requires the following conjecture. 
Namely, if a causal set C approximates to a spacetime (M,g), then (M,g) is unique, up to 
modifications to it on scales < p~ l . In other words, this conjectures that all the meaningful in- 
formation about the geometry and topology of (M,g) at scales » p" 1 is contained in the causal 
set; continuum information below these scales is irrelevant since the discrete substructure, i.e., 
the causal set, is fundamental. 

From a purely mathematical point of view, this conjectural is very intriguing. While it has 
been verified in several different cases a general proof is still not known, though considerable 
progress has been made in this direction [5J [6]. A key question is how to extract continuum 
topological and geometric properties from C using purely order theoretic information. Unique- 
ness of the approximating spacetime with respect to a given geometric or topological property 
then follows, i.e., any two spacetimes which are approximations to C must share this property 
on scales » For example, for a causal set C that is approximated by flat spacetime, the 
Myrheim- Meyer dimension gives a good estimate of the spacetime dimension [TJ [8] , while the 
length of the longest chain or total order between elements in C gives a good estimate of the 
time-like distance [9j. An estimator for spatial distance in this case has also been obtained 
|10| . Additionally, the homology of spatial hypersurfaces can be constructed from the causal 
set underlying a globally hyperbolic spacetime [111 I12j . A very important recent result is the 
construction of the scalar curvature from which the causal set action is obtained [13] . 

A natural question to ask of the flat spacetime results of [HJ [9] is how they are modified 
in the presence of curvature. While it is true in the continuum that a sufficiently small neigh- 
bourhood of a point is approximately flat, the corrections from curvature can in fact be well 
quantified using Riemann normal coordinates in a convex normal neighbourhood of any point. 
However, we run into the following issue in the discrete case: there is at present no known purely 
order theoretic definition of a "small" neighbourhood of an element in a causal set which cor- 
responds to a convex normal neighbourhood. At best it may be possible to find approximately 
flat subsets in a causal set but it is still unclear how to do this in a systematic way |14j . 

This issue however will not be the focus of our current investigation. Instead we consider 
only those A r -element causal sets C which are approximated by a small causal diamond or 
Alexandrov interval I[p, q] in a generic curved spacetime (M,g), where the smallness parameter 
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is given by the proper time T between the events p and q. We expand the metric in Riemann 
normal coordinates (RNC) about an origin r = (0, . . . , 0) 

9ab{x) = ^(O) - ^x c x d R acbd (0) + 0(x 3 ), (3) 

where the 0(x) correction to the metric vanishes since r^ 6 (0) = 0. The RNC has been used to 
calculate the volume of a small interval I[p, q] [7| 115} [T6] for which the first order correction to 
the volume of I[p, q] due to the effect of curvature occurs at 0(T 2+n ). 

The starting point of our analysis is Meyer's work [8] in which a general expression was found 
for the average abundance of ^-chains or /c-element total orders in a causal set Co which is 
approximated by an Alexandroff interval Io [p,q] in flat spacetime. Following Myrheim, Meyer 
used this to find a dimension estimator for the dimension of Io\p, q] employing only (Ci), the 
average abundance of elements in Io[p, q] and (C2}, the average abundance of relations. Using the 
RNC expansion and with the help of [16] we extend this analysis to the curved spacetime case. 
We find that the in I\p, q] satisfy a recursion relation and depend on the scalar curvature 
-R(O) and the time-time component of the Ricci tensor Rqq(0). It is then an easy exercise to 
invert these relations and find expressions for R(0) and -Roo(O) in terms of (Ci), (C2), (C3) for 
fixed spacetime dimension. We find that the Myrheim-Meyer dimension estimator is insufficient 
to determine dimension in a generic curved spacetime and construct a new dimension estimator 
for generic non-flat spacetimes using (C±), (C2), (C3), (C4). 

As mentioned earlier, the scalar curvature of an element in a causal set was first calculated by 
Benincasa and Dowker |13] using a curved spacetime expression for a non-local D'Alembertian 
on a causal set. They showed that R(q) for an element q in the C can be expressed in terms 
of the abundance of k-element "inclusive-intervals" which are order theoretically very distinct 
from fe-chains. Like fc-chains they too have a bottom element e\ and a top element e^, but 
unlike /c-chains, every element in the interval I[ei,efc] which is the intersection of the future of 
ei and the past of belongs to the inclusive interval, and has precisely k — 2 elements satisfying 
e i -4 &i ~< Cfc. In contrast, the expression for R(0) that we find depends only on the abundance of 
fc-chains. This may suggest that for manifold-like causal sets there are hidden relations between 
these seemingly different order theoretic entities. 

The plan of our paper is as follows. In Section [2] after first presenting some basic definitions, 
we reproduce Meyer's results for the (C^) thus setting the notation that we will use in the rest 
of the paper. In Section [3] we present the main calculation in the paper, where we use the 
RNC to obtain the lowest order curvature correction to (C^). We show that there is a recursion 
relation between the coefficients in the expression for the different (Cfc), but that the form of 
the dependence on the R(0) and i?oo(0) is the same for all k. In Section |3] for a fixed spacetime 
dimension we find expressions for R(0) and i?oo(0) which depend only on (C\), (C2), (C3). The 
coefficients in these expressions again have a simple dependence on dimension. In Section [5] we 
point out that the Myrheim-Meyer dimension estimator is insufficient in a generic curved space- 
time and find a new dimension estimator using (Ci), (C2), {C3), (C4} . An important question 
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in these calculations is how the error decreases with the sprinkling density p (i.e., the inverse of 
the volume cut-off) - the larger p is the closer one comes to the continuum. In Section [6] using 

2k — 1 

the technique developed in [8] we show that the error in (Cfc) grows as p~^~ which means that 
the error in R(0), i?oo(0) and n goes like p _1//2 , thus going to zero in the continuum limit. We 
discuss the implications of our results in Section [7] and the questions that need to be addressed 
in the future. Finally, the Appendix contains explicit calculations of the results that appear in 
the main body of the paper. 



2 The Abundance of /c-Chains in Flat Spacetime 

A /c-chain in a causal set C is a fc-element total order, i.e., a set of elements {e±,e2, ■ ■ ■ , e^}, 
e.j E C such that e% -< e^+i for all i. For any finite element C, the number Cf. of fc-chains 
is therefore invariant of the choice of labeling of C. This makes Ck a good observable. Note 
that the and ej+i need not be "linked", i.e., there could exist an element e G C such that 
ei < e < Moreover, e± ~< e ~< does not imply that e belongs to the fc-chain. In contrast, a 
fe-inclusive interval is defined as ifc[ei,efc] = Future(ei) n Past(efc) [15] . Along with the elements 
ei, efc it also contains precisely /c-elements. However, every e € C such that e\ ~< e ~< belongs 
to /jfc[ei,ejfc], which means that the order theoretic structure of a A;-chain is very different from 
that of a fe-inclusive interval. The fact that one can express the discrete scalar curvature both 
in terms of the abundance of the inclusive intervals as shown by Benincasa and Dowker [13] 
and in terms of the abundance of /c-chains as we will show in Section U] thus suggests a hidden 
connection between the two. 

We now reproduce Meyer's results in n-dimensional flat spacetime (Mo, 77) using notation 
that we will find convenient in the curved spacetime generalisation. Let p,q E Mo such that 
p = (— T/2, 0, . . . , 0) and q = (T/2, 0, . . . , 0). For a causal set Co that is approximated by Io[p, q] 
for a given sprinkling density p, the average abundance of elements, or 1-chains {Ci) v is given 
by 



(C±) v = pV = p j dxi 



lo [p,q] 

2p I dt I dr r n ~ 2 [ dn n _ 2 = p 2An ~ 2 (lY = p( T n , (4) 





?n-2 



n(n - 1) V 2 



where A n _2 is the volume of the unit (n — 2) sphere S n ~ . Next, the average number of 2-chains 
or relations in Io\p, q] is given by the probability of there being a pair of elements x±, x-i E lo\p, q] 
such that x\ -< X2, or 

(C2) v = P 2 dxi \ dx 2 (5) 



Mp,q] J+{xi)nI [p,q] 
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Recognising that the integral over dx 2 is simply the volume of the smaller interval Io[^i, q] and 
using EqfllJ 



2 ^n-2 f , nn _^ a r(n + i)r(§) 



x ln H 2 n n(n-l) J 4T(^) 
Io [p,g] 

Meyer was able to similarly use the nested integral expression for {Ck)n 

{C k ) v = p h J dx x J dx 2 ... J dx k = p k J dx 1 (C k - 1 (x 1 )) v (7) 
iofp,<?] J+(xi)ni [p.g] J + (xk-i)rii \p,q] iob,g] 
to find by induction the general form 

(Ck) v = P h XkV k = p k C k T kn , (8) 

where 

lf Tin + l) ^- 1 T(§)r(n) _ / 2A n _ 2 \ k k 

Xk= k{-^—) r( M )r ((^)' Cfc= V2«n(n-l)J = Co ^ 

with defined as in Eqn (j4|). Note in particular that xi = 1. We will find it useful to express 
(Ck}ri as 

(C k ) v = P k Ck-i J dx x T[ k - 1)n = - l)n) (10) 

lo[p,g] 

where Ii(m) is evaluated in Eqn (1T5f) of the Appendix. As discussed above, the average number 
of chains in a finite element causal set C is itself a covariant observable. In particular, the 
distribution of the abundance of /c-chains as a function of in a finite element causal set can 
be compared with the distribution of (Cfc)„; if the two distributions agree, it is an indication 
that the C may be well approximated by flat spacetime and is therefore manifold-like. A similar 
comparison using /c-inclusive intervals was found to be useful in determining flat spacetime 
behaviour in a model of 2d causal set quantum gravity [17j . It is therefore important to find a 
generalisation of (Ck) v to curved spacetime. 

Meyer obtained a dimension estimator from (C^)^ by observing that the ratio 

f (n) _ (c 2 ) v _ r(n + i)r(3) 
fo(n) = ~m ~ 4r ( f) 

is only a function of n. Thus, one has an expression for the dimension which depends only on 
order-theoretic information in the causal set. Indeed, fo(n) is one-half of Myrheim's ordering 
fraction 
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where R = (C2) is the number of relations and N = (Ci). In 2 spacetime dimensions for 
example, f(2) = 1/2, i.e., the inverse of the spacetime dimension. In Section[5]we will show that 
Eqn (jlip does not suffice in curved spacetime and there is need for a new dimension estimator. 



3 The Abundance of k- Chains in a Small Causal Diamond in 
Curved Spacetime 

The RNC expansion to order T 2 gives an expression for (C\) in I[p, q] [7] 115j 



(d) = P V = pJ d Xl = pV (l + ai R(0)T 2 + fcRvo^T 2 ) 

l\p,q] 



where 



OLl = — ■ 



24(n + l)(ra + 2) : 
and which uses the RNC expansion 



ft 



n 



24(n + 1) 



V=gl = 1 - -x^R^O) + 0(x 3 ). 



Now, the average number of 2-chains or relations is given by the similar generalisation 
(C 2 ) = p 2 J y/-gi dxi J yj-gi dx 2 = p 2 j V~9i dx{Vi, 



(13) 



(14) 



(15) 



(16) 



J+(xi)r\I\p,q\ 



where Vi denotes the volume of the region J + (xi) n I[p, q]. Using the covariant form of Eqn 
(fl~3l) we see that 

<C 2 ) = p 2 Ci J dx!l?(l- ^xtxlR^O) + a 1 T 2 R(y 1 ) + ^T^R^(y 1 )^+0(T n+3 ) 



p 2 Ci 



+ J dxxTff—xIxIR^O) + ai T 2 R( yi ) + PtTfTfR^yt)) 



0(T 



n+3\ 



lob, 9] 



(17) 



where we have split the integral in the manner of [7J : the first is a flat spacetime integral over 
the curved spacetime interval I[p, q], whereas the second is the contribution from the curvature 
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p 



Figure 1: An Alexandroff interval I[p, q] in flat spacctimc. T\ is the proper time between the events x\ 
and q and yi is the midpoint of 5]. 

terms over the flat spacetime interval lob, l]- Vi ls the midpoint of the interval l[xi, q] as shown 
in Fig [TJ Using the light-cone coordinates u = t — r , v = t + r we see from Fig [1] that 



if -(I)*-* *i-v(f -*)(!- 4 (18> 

Thus, 

T?T?R^( yi ) = ^i?oo(yi) - x^TRo^) + x^R^iVi)- (19) 
The first integral was evaluated in [16] and shown to be of the form 



dxiT?= J dxxT^fl + ^RooiO)) 



(20) 



i[p,g] io[p,<?] 



for any non-negative integer m. Moreover, as can be readily seen, to order T 2 , R liv {y\) can be 
replaced with R^ u (0). We can moreover simplify the expressions in (C2) substantially by using 
the symmetries of Iq [p, q] ■ Expanding the term 



/ 



dx 1 x^x u 1 R IJtv {Q)T^ = j dx 1 t 2 1 R 00 (0)T[ n +2 J dx 1 tx\R Qi (0)TJ n + J dxxx^Ri^T^ 



iob,<?] lob-?] io[p,<?] iob>9] 

(21) 
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for m a positive integer, we see that the cross-terms do not contribute, so that we are left with 
dx^x^R^T? = i? 00 (0) / dxitlT^ + Y^RiM / d Xl (x{) 2 T^, (22) 

i t r i 1=1 t r i 



lo [p,q] 



lo \p,q] 



iob,g] 



since the last integral is independent of the spatial direction i, due to the symmetry of lob, 
Gathering the coefficients of R(0) and -Rqo(O) 



<c 2 ) = P 2 ci 



^oo (0) 



h(n) + R(0) 

-rp2 



{Pi - hh(n)+a 1 li(n + 2) 
b 



+ 



>A h(n) + (A - i)I 2 (n) - Il 3 (n) + /3 1 I 4 (n) 



(23) 



n-l 



where we have used Yl ^ii(O) = -Roo(O) + -R(O) and we define the general class of integrals 

i=l 



dx\T] 



lo [p,?] 



1 ■ 



I 2 (m) 



dxiT^rl cos 2 6»i 



lo [p,<?] 



I 3 (m)= / dx x t{Tl\ I 4 (m) 

lo [p,<?] 



(24) 



lo [p,q] 



for non-negative integers m. These integrals have been evaluated in the Appendix. Using the 
terminology defined therein, we will find it useful to re-express the last three integrals in terms 
of the first 



and 



I 2 H = / 2 (m)I 1 (m)T 2 , I 3 (m) = / 3 (m)I 1 (m)T 2 , I 4 (m) = / 4 (m)Ii(m)r 2 



l 1 (m + 2)=g 1 (m)I 1 (m)T 2 , 



where #i(m), / 2 (m), /3(m) and / 4 (m) are defined in Eqn ([ST]) of the Appendix . Given Eqn. 
(fTUl) we find that 



(C 2 > = <C 2 >, 



1 + T 2 a 2 R(0) + T 2 /? 2 i? 00 (0) 



(25) 



where 



«2 = (A - ^)h{n) + aigi(n) 

D 

ft = ^ + (/3i-^)/2(n)-i/ 3 (n) + /3 1 / 4 (n). 



(26) 



S 



We find that 



«2 



in 



24(2n + 2)(3n + 2) 
in 



24(3n + 2) ' 
We can go one step further and calculate 

(C 3 > = p [ dxi yf^i(C 2 { Xl )) 



P% 



h(2n) + R(0) ((fa - i)I 2 (2n) + a 2 I 1 (2?i + 2)) + 



T 1 



^oo(O) ( ^Ii(2n) + (/3 2 - -)I 2 (2n) - -I 3 (2n) + /3 2 I 4 (2n) 



where again we have used Eqn (110(1 and 



l + T 2 a 3j R(0) + r 2 /33^oo(0) 



Q!3 
ft 



(ft - -)/ 2 (2n) + a 25 i(2n) 



6ri 



24(3n + 2)(4n + 2) 



111 6r7 

_ + (/3 2 _ -)/ 2 (2n) - -M2n) + fafA2n) = r. 

24 KH 6 JJ y ' 6 J v ^^ 4 V > 24(4n + 2) 



This suggests an iterative formula 



with 



Wk ~ -)/ 2 (A:n) + a k gi(kn) 
b 



ft 



12(fcra + 2)((fc + l)ra + 2) 
nA; 



12((fc + l)n + 2) 



Lemma 1. To the lowest order correction in the flat spacetime expression, the average 
of k-element chains in a small causal diamond is 



(Ck) = (Ck) v 



l + T 2 a k R(0)+T 2 f3 k R 00 (0) 



+ 0(T 



where a k and ft are given by Eqn (E7 



Proof: We will prove this inductively. We have already shown it for k = 2. Now, let us assume 
(Cfc) is of the form Eqn. (|32p . Just as in the flat spacetime case, one has nested integrals so 
that 



(Cfc+i) = p j dxi^/^g{(C k (xi)) 

l\p,q] 

= P k+1 (kJd Xl T? n (l + ^fioo(O) - ixfx^(O) + a k T?R(0) + /3 fc Tf T^R^OU 



lob.?] 
+ 0{T kn+3 ) 



(33) 



where we have used Eqn ([20]) to reduce the integral over I\p,q] to one over Io[p, g] to order 
0(T kn+2 ). Using the integrals Ii,2,3,4( 7Ti ) from the Appendix, we can reduce this to 



(C, 



fc+i 



p k+1 Ck h(kn) + R(0)^ k - ^)h(kn) + a k h(kn + 2)j + 
( ^li(fcn) + (A - i)Ia(fcn) - jU(fcn) + ^(H 



1 + T^(0)( (/3 fc -i)/ 2 (fcn) + a fc5l (A ; n + 2) ) 



G 



1 



24 



Writing it in the form Eq (|32p . with 



T 2 R 00 (0) — + (& - -)/ 2 (fcn) - -/ 3 (fcn) + /9 fc / 4 (*n) 



6 



6' 



(34) 



a fc+ i = (J3 k - -)h(kn) + a k g 1 (kn + 2) 

O 

Am = ^; + (Pk-bf2(kn)-]:Mkn) + f3 k U(kn). 
24 b b 



(35) 



we find the desired form for (C k+ i). It then follows from the expressions for gi(m), f 2(171), f 3(171) 
and / 4 (m) (Eqns (fTH]) : (f75]) , ([7^]) . t[5H|) ) that 



n(fc + 1) 



12((k + l)n + 2)((k + 2)n + 2) 
n(k + 1) 



12((fc + 2)n + 2) 



(36) 
□ 
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4 Scalar Curvature from the Abundance of k- Chains 



For a fixed n (C k ) contains three unknowns, T,R(0) and Rqq(0). Thus, we need at least three 
values of k in order to determine R(0). For each {C k ) the lowest order correction due to curvature 
is 0(T kn+2 ). Hence, as in the flat spacetime calculation of the Myrheim- Meyer dimension Eqn 
(fTT|) we must take appropriate powers of (Ci), (C 2 ) and (C3) to be able to compare their lowest 
order corrections. Defining 

3 / k 1 / in, \ \ 3 / k 

Qk 



for k = 1,2,3 



Qi = T in \l + 3 ai R(0)T 2 + 3(3iR 00 {0)T z j + 0{T 6n+z ) (38) 
Q 2 = T 3n (l + ^a 2 R(0)T 2 + ^ 2 R 00 {0)T 2 ^j + 0(T 3n+2 ) (39) 
Q 3 = T 3n (l + a 3 R(0)T 2 + (3 3 R 00 (0)T 2 ) + 0{T 3n+2 ). (40) 



Thus, the Q k are independent of the sprinkling density p and hence can be used to construct 
continuum geometric parameters. It is useful to gather a few identities and definitions before 
we proceed: 

hotl x = -(kn + 2) 

*fc = Oikfik+i ~ Oik+iPk = -na k a k+ i 
k 

k + 1 

K k = ((k + l)n + 2)Q k , 

Jk = (kn + 2)K k (41) 

We eliminate the Rqq(0) term from Eqn (|38[) and (|39p and subsequently from Eqn (|39p and 
(|40l) to get the pair of equations 

o /-> \m~3n _ ^ , 3 ,T, .rr2- 



^Qi-PiQ2)T- 6n = ^ 1 + -^ 1 T 2 R(0) (42) 
^j-Qi- ^Qsjr- 3 " = $ 2 + * 2 T 2 R(P). (43) 
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Since both R(0) and T are unknowns, we first eliminate the R(0)T 2 term: 



y* 2 Qi - 3(01*2 - /3 3 *i)Q 2 + fo^iQs )T~ in = -* 2 *1 " *1$2- 
We find after some algebraic manipulation that 

r3n = ^^i-2J2 + J 3 ). 



2n 2 

with the Jj's given by Eqn. (|4ip . We thus obtain the expression for the scalar curvature 



p/n'v 2(n + 2)(2n + 2)(3n + 2) ^ 

i?(0) = re 3 T 3n+2 (^1 " 2 ^2 + ^3) (46) 



or more explicitly 



R(0) = -2(n + 2)(2n + 2)(3n + 2)2 3 &n&- 1 {Kl ~ 2 * 2 + K *2 +a . (47) 

(J!-2J 2 + J 3 )~ 

We may additionally solve for i?oo(0) by eliminating R(0) from Eqn (f38l) . (j39j) (f40j) : 

YQi-aiQ 2 )T- 3 ™ = 0! - I^T 2 ^) (48) 
^Q2 ~ a 2 Q 3 )T~ 3n = ®2-*2T 2 R(0). (49) 



Either equation along with Eqn (j45J) gives 

^oo(O) = - 4(2n n + 3 ^ 3 + n 2 + 2) ((n + 2)Qi - (5n + 4)Q 2 + (4n + 2)Q 3 ). (50) 
As a check, let us consider the case -Roo(O) = so that from Eqn (j38|) - (|4*0|) we see that 

R(0) > =>- Q x < Q 2 < Q 3 
R(0) < Qi > Q 2 ><5s 

12(0) = Qi = Q 2 = Qs = T 3n . (51) 



Moreover, 



Kt - 2K 2 +K 3 = ; -— ^ -i?(0)T 3n+2 

d 2(ra + 2)(2n + 2)(3n + 2) y J V ' 

which therefore has the opposite sign to R(0) and is zero when R(0) = 0. 

As a further check, we note that if both i? o(0) = and R(0) = 0, Q\ = Q 2 = Q3 = T 3n 
so that not only is K\ - 2K 2 + K 3 = (2n + 2)Qi - 2(3n + 2)Q 2 + {An + 2)Q 3 = 0, but also 
(n + 2)Qi - (5n + 4)Q 2 + (4n + 2)Q 3 = which appears in Eqn {5QJ). 
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5 A New Dimension Estimator for Curved Spacetime 

As one can guess by now, the ordering fraction or equivalently the function f(rt) in curved 
spacetime clearly involves the curvature contribution non-trivially. Expanding to order T 2 the 
curved spacetime version of the Myrheim- Meyer dimension estimator is 

f(n) = ^ 

= f„(n)( 1 + T 2 (a 2 - 2 ai )fl(0) + r 2 (ft - 2ft)flco(0) ) + 0(T 3 ) 



= f ° ( "' (' + !2(n + "ikL + 2) ((dh) - fiM(0> )) + ° (r3) - <53> 

In the special case that R(0) = .Roo(O) = 0, f(n) « f (n) up to order T 2 . For a generic spacetime 
f(n) it is however insufficient as a dimension estimator and we must find a replacement. Given 
that along with n there are 4 unknowns to be solved in terms of the (Cf~), the simplest way to 
do so is to include (C4) in our analysis. 
We define 

f v 4/fc 

s fc = (54) 



analogous to the in the previous section, with k = 1, 2, 3, 4 which again is independent of the 
sprinkling density. To the lowest order correction we have the four equations 

Si = r 4n ( l + 4aiE(0)T 2 + 4£i.Roo(0)T 2 



5 2 = T 4n 1^1 + 2a 2 R(0)T 2 + 2f3 2 R 00 (0)T 2 

5 3 = T 4n (l + ^a 3 R(0)T 2 + ^ 3 R 00 (0)T 2 

S A = T An (l + ai R(Q)T 2 + (3 4 R 00 (0)T 2 ^ (55) 
Eliminating R(0)T 2 from the above we get 

= T in (e 1 -2^ 1 R 00 T 2 





- a\S 2 


7^a 3 S 2 


- a 2 S 3 


-a 4 S 3 





e 2 - U> 2 R m T 2 



T m (e 3 -* 3 R 00 T 2 \ (56) 
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from which we may eliminate Rqq(0)T to get 



3 (ja 2 Si-aiS 2 J* 2 -2^-a 3 S 2 -a 2 S3j*i = r*M -9itt 2 - 2e 2 *i 

I 0a 3 ^2 - a 2 S 3 ) * 3 - ^a 4 5 3 - a 3 S 4 ) ¥ 2 = T 4n Qe 2 * 3 - 9 3 ^ 2 ) • 

This gives us an expression for T 4n (which we check reduces to Eqn (j45j) ). After some algebra 
this gives the following implicit form for the dimension: 

(n + 2)(2n + 2)5i - 3(2n + 2)(3n + 2)5 2 + 3(3n + 2)(4n + 2)S 3 - (4n + 2)(5n + 2)5 4 = 0. (57) 

Importantly, in the absence of curvature the Sk are all equal and the left-hand side is identically 
zero. Hence this cannot be a replacement for the Myrheim-Meyer dimension in flat spacetime. 
Using Uk = (k + 2)((k + l)n + 2)Sk we may write the expression more succinctly as 

U x - 3U 2 + 3U 3 -U 4 = 0. (58) 

It is interesting to note the appearance of the binomial coefficients (— l) fe ( r 7 ) for r = 4 in the 
above expression for the dimension estimator, as well as in the expressions for R(0) and T for 
r = 3. 

Since the Sk themselves explicitly contain dimension information via the Cfc it is more useful 
to expand the expression to 

(n + 2)(2n + 2) ( ^] - 3(2n + 2)(3n + 2/ ^ ' 



or defining 



Xl J \ X2 

+3(3n + 2)(4n + 2)(^^J - (4n + 2)(5n + 2) (^-) = (59) 

_ . fc _i/ 3 \ (fcn + 2)((fc + l)n + 2) 

Wfc = H U-v ( } 

we get our final expression for the dimension estimator in curved spacetime 

4 

^u;,(n)<C fc ) 4A = 0. (61) 
k=i 



6 Calculating the Errors 

We expect that as the sprinkling density p increases, our curvature and dimension estimators 
should do a better job of reproducing the continuum results. While the geometric parameters 
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themselves do not depend on p, it is clear that the error will. The (rms) error 5Ck = y/ACj~ = 
yj (C|) — (Cfc) 2 , where AC^ is the variance. We follow the analysis in [8] to find the dependence 
of ACfc on p in the RNC. Unlike the flat spacetime case however, it is not the continuum volume 
that must be increased to improve accuracy since this region must still be "small" for the RNC 
to be valid. Instead, it is the sprinkling density p that should be increased for reducing the 
error. 

Let us begin with k = 2, so that AC2 = (Cf ) — {C2) 2 ■ Now (Cf ) is the probability of finding 
two sets of 2-chains in I\p, q], with the possibility that some of the elements can coincide. Let 
us call the points of these two chains x, x' and y, y', with x -< x' and y -< y' . Thus, (Cf ) gets 
contributions from each type of coincidence. The first is simply that there are no coincidences, 
i.e., that all 4 points x,x',y,y' are distinct, which gives a contribution (C2) 2 . Although this 
term ~ p 4 , it cancels out in the expression for the variance and therefore plays no role. The 
next type is the one-coincidence case. For this there are two types: (i) x = y' or y = x' , so 
that the two 2-chains collapse to a single 3-chain, and (ii) x = y or x' = y', corresponding to 
the probability for a three element "V" or "A" shaped causal set. The contribution from (i) is 
clearly twice that of (C3) whose p dependence is ~ p 3 , while the contribution from (ii) is 

2p 3 Jdxy / -g(x) Jdx'^-g(x') J dy'^-g{y') = 2p j dx{d{x)) 2 ~ p\ (62) 

I[p,q] J+{x)nI\p,q] J+(x)nl\p,q] I]p,q] 

Finally, there is a contribution from two-coincidences x = y,x' = y', which is just (C 2 ) which 
goes as ~ p 2 . Thus, AC2 ~ p 3 or 5C2 ~ p 3 , with the dominant contribution coming from the 
one-coincidence case. 

In order to calculate the error for all our geometric parameters, we need to perform a similar 
analysis for k = 3, 4. In each such case, the dominant contribution to the error comes from the 
one-coincidence case, since the no-coincidence contributions simply cancel out. 

For k = 3, if x -< x' -< x" and y -< y' -< y", the no-coincidence case is again simply (Cf) and 
again cancels out. The one-coincidence terms include: (i) x = y" or y = x" which is the 5-chain, 
(C5) ~ p 5 , (ii) x = y' or y = x' or x" = y' or y" = x' which contribute 



V Jdyyf^) J dxy/^gJxjiCiix))^)) ~ p 5 , (63) 

l\p,q] j+(y)nl\p,q] 

(iii) x = y or x" = y" which contribute 



2p dx^f^){C 2 {x)) 2 ~ p\ (64) 



and finally (iv) x' = y' which contribute 



p 3 dxj^gjx) dy^g~{y) / dx' ^(^{dix)) 2 ~ p\ (65) 

Jl\p,q] Jl\p,q\ J 

J+(x)nJ+(y)nI\p,q] 
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Figure [2] shows the contributions to (C|) from the case of no-coincidence and the case of 1- 
coincidence. Since the no-coincidence term cancels out, 5C^ ~ p2. 




Figure 2: The contributions to (Cf ) from no-coincidcnces and one-coincidence are shown. 

A similar analysis shows that for (Cj) the no-coincidence term cancels out and the one- 
coincidence cases lead to a dependence ~ p 7 , so that SC^ ~ pa This generalises in a straightfor- 
ward manner to 

bC k ~p^. (66) 
We are now in a position to calculate the dependence on p of the errors in the Q k and S^: 

SQk ~ ^3 {C k f' k - l 8C k ~ p" 1 / 2 . (67) 

<5S fc ~ J-^V*-!^ „ p -i/2 (68) 

This immediately means that the errors 5T 3n , 5R(0), 5Rqo(0) go as p" 1 / 2 and hence become 
smaller as p increases. The error in the dimension estimator is similarly given by 



and 



8nl (4n + 6)Si - 3(12n + 10)5 2 + 3(24n + 14)53 - (40n + 18)5 4 

(n + 2)(2n + 2)5S 1 - 3(2n + 2)(3n + 2)552 + 

3(3n + 2)(4n + 2)<55 3 - (4n + 2)(5n + 2)<5S 4 ^) 
=> <5n ~ p~ 1/2 . (69) 

7 Conclusions and Remarks 

In this work we have found expressions for the proper time T (Eqn (|45p ). the scalar curvature 
.R(O) (Eqn (07])), the time-time component of the Ricci tensor -Roo(O) (Eqn ([50]) ). and a new 
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dimension estimator (Eqn (|6ip ) from a causal set underlying a small causal diamond I[p, q] in a 
generic spacetime in arbitrary dimensions. We find that the errors in these estimators goes as 
/9 -1 / 2 , thus becoming smaller as the sprinkling density is increased, while keeping the volume 
of q] fixed. Our results not only verify the deep relationship between order and Lorentzian 
geometry, but also provide new observables that can be used to assess whether a causal set is 
manifold-like or not. 

Our calculation moreover brings to light an intriguing connection between two seemingly 
disparate order theoretic structures in a causal set. While our expression for the scalar curvature 
is purely in terms of the abundance of /c-chains, the Benincasa-Dowker(BD) scalar curvature Rbd 
|13j is constructed from the abundance of /^-inclusive intervals. In 4-dimensions for example, 
their expression for the scalar curvature is 

Rbd(0) = -^(1 - (N 2 (0) - 9iV 3 (0) + 16JV 4 (0) - 8JV 5 (0))) (70) 

where iVfc(O) are the number of k- inclusive intervals Ik[x, 0], where x -< (@. As noted in Section 
[2j /c-chains and fe-inclusive intervals are very distinct order theoretic structures. The continuum 
geometry however seems to link them via the scalar curvature. If such a relationship exists, 
does it for example indicate manifold-likeness in a causal set? This and several related questions 
remain to be investigated. 

Further comparisons between the two expressions for R(0) are warranted. Both contain 
alternating sums, although the coefficients differ markedly. For one, the BD expression appears 
to have a strong dimension dependence in the number of terms required - for n = 2 the sum 
is truncated at k = 4, while for n = 4 it is truncated at k = 5. Our expression for -R(O) is 
in this sense independent of n - it requires (C\), (C2), (C3) in all dimensions^]. Moreover, the 
systematic determination of the coefficients for Rbd for arbitrary n is fairly involved, whereas 
the coefficients in Eqn (j47j) and Eqn (I58P are simply the binomial coefficients (— and 
(— l) fc (^), respectively, in all dimensions. 

Apart from these, there are deeper differences in the two expressions for R(0). Rbd(0) is 
constructed from all inclusive intervals of the form Ik[x, 0] in the causal set for k = 1, . . . 5. It is 
therefore an essentially non-local expression since it depends on the structure of the causal set 
throughout a past (or future) neighbourhood of the element in the causal set. The expression 
Eqn (|47p on the other hand is only valid in a small causal diamond, and hence is strictly local 
and dependent on a proper choice of neighbourhood of the element in the causal set. This means 
that the BD form for -R(O), being "neighbourhood independent", can be readily used to obtain 
an action for the entire causal set, which has strong implications for causal set theory. Unless 

2 Our notation differs from [13] where k is replaced by k — 1, in keeping with our definition of a fc-chain. 

3 As may be already evident to the astute reader, one could replace these three values k = 1, 2, 3 with any other 
fci,fe2,fe3, to get an expression for -R(O) in terms of (C^), (Ck 2 ), (C* 3 )- What is important is that this choice is 
not dimension dependent - every choice works for every n. 
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the definition of a small neighbourhood can be made entirely order theoretically in the causal 
set, our expression for R(0) on the other hand cannot be used to obtain an action in a simple 
manner. Of course, in the specific case when the entire causal set is approximated by a small 
causal diamond q] the action is simply 



where e s denotes an element in C, since to this approximation R(x) is the same throughout 



Nevertheless, the geometric estimators we have found in this work take a big step towards 
the order-Lorentzian geometry correspondence. We have a new dimension estimator which 
can determine the manifold dimension for a curved spacetime with greater accuracy than the 
currently available flat spacetime Myrheim- Meyer estimator. The local nature of the expression 
for the scalar curvature can also be seen as an advantage since only a small neighbourhood of 
an element in the causal set is required to determine the curvature, rather the than its entire 
past. An obvious next step is to follow up with a numerical analysis of causal sets which are 
approximated by different curved spacetimes and to see how well our estimators work in these 
cases [18] . 

These observables can also be used as additional tests of manifold-likeness of a causal set. 
In [17] the expectation values of a range of such observables, including the abundance of 
/c-element inclusive intervals was found from Monte Carlo simulations of 2-d causal set quantum 
gravity. Comparing with the flat spacetime distribution of the Nk, these observables were used 
to demonstrate that the dominant contribution to the causal set path integral in 2d comes from 
flat spacetime suggesting that manifold- like behaviour is emergent. The estimators we have 
obtained in this current work could well be employed to calculate more accurately how close to 
flatness one is - is the 2-d universe just a little positively curved, for example? 

The extended hope of the current analysis is also that as more of the order theoretic basis 
of geometry is uncovered it may be possible to find a meaningful order theoretic definition of 
locality which translates to our commonly held (albeit Riemannian geometry based) notions of 
locality in the continuum. 



N 




(71) 



s=l 
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8 Appendix 



To evaluate the (Ck) we find it convenient to transform from the Cartesian coordinates x % to 

„„i, ■ „i i„„ j- „j.„„ „j ti(n\ ...1 



-LU evaiucue me we mm n cunvemeiii iu 

spherical polar coordinates x l = rf l (£l), where 

f(e 1 ,9 2 ...9 n _ 2 ) = < 



( n—i—l 

n sin 8 k cos 6 n -i (i > 1) 

k=l 



n-2 

n sin6i fc 

fc=l 



(i = l) 



It is also useful to express the radial coordinate for any aii G I[p, q] in light-cone coordinates: 



«i - u\ 1 / YT 
ri = — 2 — = lU^"" 1 

k=o v / v / 
Since the metric r/(u,v) in Iq[p, in light cone coordinates is 

2 ( v ~ 

ds = —dudv + I — - 

the associated measure of any integral over Io [p, <?] is 

^n— 3 



T 

Vi 

2 

'T 



(72) 



i / \ n-2 /n-3 



8=1 



.. /n— 3 \ n— 2 

=^(n- ( "- 2, ".)S(-i) v , , V2 V2 

v i=l 7 fc=0 V 7 V 

Finally, we will find it useful to define the following pair of integrals: 
L(a,v,T) = 



r 



(73) 



T 
" 2 

J 

2 



'T 



L(a. 



f / T 1 \ b T^a+b+2 

b,T) = dvl v )L(a,v,T) = - - ■ -. 

' J \2 J K ' ' ' (b + l)(a + b + 2) 

T 

2 



(74) 



2 

The above identities will now be used to evaluate the set of integrals Ii,2,3,4( r ^) required for 
the calculation of (C k ) (Eqns (fTTD . (1281) (PD )- 
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him) = J dx x T^ 

T 

2 Vi 

= J dv\ j dui J dVLyJ —r]{u\, v\) T"' 



T _T 
' 2 2 



l — n \ / 



fc=0 

An _ 2T n+m (n-2)!(f)! 

(n + m)(n + f -1)!" 

2. For m + 2 the above equation gives us the relation 

3. As noted in Section [3] any integral of the form 

h(m)= J dxMfT? 1 
io[p,g] 

is independent of the spatial direction i because of the spatial symmetry of Io\p, q]. 
we can choose x n ~ l = rcos#i to simplify our calculation so that 



h(m) = J dx{T?r\ cos 2 0i 



dvi dm dil 2n+1 \2~ Vl ) [ Vl ~ Ul 



T _T 
' 2 2 



A n _ 2 T n+m+2 n!(f)! 

2 n+i ( n _ i)( n + m + 2) ( n + f + 1)! 

T / "(^ + m ) 

h{m)l 



4(n + m + 2)(n + f + l)(n+f ) 
= him)T 2 him). 

The next integral can be split into three parts 



him) 



J d Xl t\T? = If + 1| + 1§ 
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where, using a^ = n — 2 — k + ^ ,bk = k + ^ 
li{m) = = l -fd Xl u\T™ 

k=o V / V / 

(80) 



4(m) = ^Jdx.vfT™ 



n-2 



2 n+l 

fe=0 



^(-l) fc ( n fc 2 ) (^L(a fc , b fc> T) - TL(a fc , b fe + 1, T) + L(a fc , b + 2, T)) , 

(81) 



and 



mm) = \f d XlUlVl Tr 

t — n \ / \ 



A ™" 2 
At- 2 

2^ 

fc=0 



~L(a fc ,b fc + l,r)+L(a fc + l,b + l,r) 



After some algebra, we find that 



j / \ t / \t , 2 8n + m(m + 2)(m + n + 2) , t2 ,, v , a „* 

I 3 (m) = Ii(m)T — — u — 777-— — rr = Ii (m )T / 3 (m). (83) 

4(2 + m + n)[2n + m)(2n + m + 2) 



4. Finally, 



I 4 (m) = / cfci( --ti ) Tf 



T 



rp2 

—h(m)+h(m)-Th(m). (84) 



Evaluating 

I 4 (m) = J dxihT r ' 



m + n (f + l)(m + ra) 



2 V m + n+1 (m + n + 1)(^ + n) 
h(rn)TU(rn). (85) 
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We can then express 



U(m) = li(m)T 2 f 4 (m) 



(86) 



where f 4(171) = 3 + U(m) + h{m). 

Gathering these expressions 

(ra + m)(f + 1) 



/2H 



n(n + m) 



(n + m + 2)(n + ^) 



4(n + m + 2)(n + f + l)(n + f ) 
^ - / 4 (m) + h{m). (87) 



/3W 



8n + m(m + 2)(m + n + 2) 



/4H 



4(2 + m + n)(2n + m)(2n + m + 2) 
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